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Abstract 

We study the non-equilibrium dynamics (purely dissipative and relaxational) 
in a semi-infinite system following a quench from the high temperature dis- 
ordered phase to its critical temperature. We show that the local autocorre- 
lation near the surface of a semi-infinite system decays algebraically in time 
with a new exponent which is different from the bulk. We calculate this 
new non-equilibrium surface exponent in several cases, both analytically and 
numerically. 
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There has been a lot of current interest in understanding the growth of correlations in 
a system after being quenched from the high temperature disordered phase to or below 
its critical temperature (T c ) [|TJ. In either case the system exhibits dynamic scaling at 
late stage of the growth. The growth is characterized by a single time dependent length 
scale. For quench to below T c , this length scale characterizes the linear size of the growing 
domains of competing broken symmetry phases. On the other hand, for quench to T c , 
it characterizes the length scale over which equilibrium critical properties are established. 
A lot of theoretical and experimental efforts have been directed in determining the time 
dependence of this length scale and the scaling of the equal-time correlation functions. It 
was, however, realized later that even the two-time correlation functions have interesting 
dynamical scaling . In particular, the auto-correlation function, measuring the memory 
of the initial conditions retained by the system after time t, decays algebraically with time 
0. This has been verified by exact calculations in a few cases numerical simulations 

TM and very recently experimentally |TT[ for quench to T < T c , in a liquid crystal system 



using video-microscopy. 

In static critical phenomena, it is well known that the critical behaviour near the bound- 
ary of a semi-infinite system is drastically different from the behaviour deep inside the bulk 
rjj-|n|. It is therefore natural and important to know whether there are similar modifica- 



tions in the dynamical behaviour near the boundary. For the critical dynamics following a 
quench to T c , it has been argued that even the semi-infinite system is characterized by a 
single time dependent length scale £ (t) which has the same growth law as the infinite system 
f20| , |2l"ll . However, in this letter we demonstrate, both analytically and numerically, that the 
temporal decay of the critical auto-correlations near the boundary of a semi-infinite system 
is characterized by new exponents different from that in the bulk. 

We consider a semi-infinite 0(n) model in the space [x = (r, z)\ which extends over 
infinite space in d — 1 directions (denoted by r) and over only positive half-space in one 
direction (z > 0). The system is assumed to be translationally invariant in the (d — 1) 
directions and this invariance is broken in the z direction due to the presence of a surface 

— * 

at z — 0. The model is described by an n component order parameter field (p = [0i, . . . , 0„] 
and a coarse grained Landau- Ginzburg free energy functional with an additional surface 



contribution p6[ , 

F{$) = \J d d x[{V$f + r fi + -^f + c5(z)ft], (1) 
where the integration in Eq. (1) is over the half-space z > 0. The equilibrium properties 



of this model have been studied in detail [15, lq]. Depending upon the value of c, different 



types of surface orderings take place. There exists a special value c = c* such that, for 
c > c*, the surface orders along with the bulk at bulk T c . This parasitical transition is called 
"ordinary" transition [|l^,|n| . For c < c* and in high enough dimensions (such that a d — 1 
dimensional surface can order), the surface orders first as the temperature is lowered while 
the bulk is still disordered ("surface" transition) and then as the temperature is lowered 
further the bulk orders in presence of an ordered surface ("extraordinary" transition). The 
value c = c* is a special point where the critical exponents are different from the ordinary 
or surface transitions. This is called the "special" transition. Within mean field theory, 
c* = but becomes nonzero for d < 4 due to corrections arising from fluctuations |16| . The 
critical exponents associated with these different types of transitions are different from each 
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other and from the bulk values fll6| . For example, an exact calculation exists [|T9| for the 
2-d semi-infinite Ising model (where the surface undergoes an "ordinary" transition at the 
bulk T c |16[]) which shows that the two-point correlator between two-points on the surface 
(separated by r) decays faster as r~ l for large r than the bulk correlator that decays as 

r -l/4 

In this paper, we consider the non-conserved dynamics of the order parameter field in 
presence of a surface following a quench from the high temperature disordered (T > T c ) to the 
bulk critical point T = T c and ask whether the presence of the surface modifies the dynamics 
near the surface. Far from the surface one should recover the critical dynamics of a truly 
infinite system for which several results are known. For example, it is now well established 
jl3| , [T2|| that the bulk equal-time correlation function, G(x,t) = {4>(x' ,t)4>(x' + x,t)) exhibits 
dynamic scaling, G(x,t) ~ x^^ 2 ^ g c (x / ^(t)) where g c is a universal scaling function and 
£(£) ~ t l / z is the time-dependent correlation length. 77 and Z are the usual static and 
dynamic exponents and the () denotes an average over all possible initial conditions (corre- 
sponding to the equlibrium distribution at the initial high temperature) and over the history 
of time evolution. The bulk two-time correlation function C(x,t) = (4>(x f , 0)<ft(x' + x,t)), 
measuring the correlation with the initial condition, also exhibits dynamic scaling |13| , |T2|1 , 
C(x,t) ~ [£,(t)]~ Xc fc(%/£,(t)) where / c (0) is a constant of 0(1). The exponent A c , charac- 
terizing the decay of the bulk autocorrelation, A b {t) = (<j)(x,0)<p(x,t)) ~ [£(i)]~ Ac , is a new 
critical exponent ||T2| , |T3ll in the sense that no simple scaling relation has been found relating 
it to other static or dynamic critical exponents. For an infinite system, A c has been calcu- 
lated analytically for the 0(n) model in the limit n — > 00 and also within e-expansion where 
e = 4 — d (d = 4 being the upper critical dimension) ||13|| . For Ising model in d = 2 and 3, 
A c has been determined numerically fll2| . 

The specific dynamical quantity that we calculate explicitly in this paper for the semi- 
infinite system and show that it gets drastically modified due to the presence of the surface, 
is the decay of the auto-correlation A(z,t) = (<f)(r, z,0)<j)(r, z,t)) with time t. In the limit 
z — ► 00, we recover, as expected, the bulk results A(oa,t) ~ [£(*)] Xb where we denote the 
bulk A c by X b . However, for small z near the surface, we find that the autocorrelator decays 
as A s (t) ~ [£(^)] -Aa where A s is a new dynamical surface exponent different from A&. Also, 
the value of X s depends explicitly on the type of the surface transition. In this paper, we 
calculate A s analytically within e-expansion and in the n — > 00 limit for the "ordinary" and 
"special" transition. Also, we determine the value of A s numerically for the two-dimensional 
Ising model. 

The model-A dynamics of the order parameter is governed by the Langevin equation, 

where F is given by Eq. [1] and ff(x, t) is a Gaussian noise with zero average and a correlator 
(r)i(x,t)r)j(x' ,t')) = 2kBTSij5(x — x')5(t — f), where T is the temperature. We first consider 
the Gaussian theory where one neglects the interaction (set u = in Eq. (1)) and which 
is valid for d > 4. We define the Fourier transform, G(k,z,z',t) = J <i d-1 (r — r')G(r — 
r 1 , z, z' ,t) exp[ik.(r — f f )] where k is a (d — 1) dimensional vector in the reciprocal space. 
Then from Eq. (2), at the critical point (r = and setting ksT c = 1), G(k, z, z',t) evolves 
as, 
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d t G(k, z, z', t) = [-2k 2 + d 2 + d z , 2 ]G(k, z, z', t) + 25(z - z') 



(3) 



with the boundary condition, d z G = cG at z = and the initial condition, G(k, z, z', 0) = 
A5(z — z') (this "white noise" form of the initial condition corresponds to quench from the 
infinite temperature where the field 4>(x, t) is completely random and A controls the size of 
initial onsite fluctuations in 0). Similarly, the symmetrized two-time correlation function 
C s (k, z, z', t) defined as the Fourier transform of C s (f, z, z', t) = l(((/)(f, z r , 0)(f>(f + f, z, t) + 
0(7^, z, 0)(f){f + r, z', t)) evolves as 

d t C s (k, z, z', t) = X -\-2k 2 + d 2 + d z , 2 ]C s (k, z, z>, t) (4) 

with the same boundary and initial conditions. By choosing the basis function 

i/>(u, z) = -^=[exp(iuz) - ° c + ™ u exp(-iuz)], (5) 

it is easy to see that the solutions to Eq. (3) and (4) are given by, G(k, z, z' ,t) = 
J^oo duip(u, z)ip* (u, z') fi(k, u,t) and C s (k,z,z',t) = dml>(u, z)ip* (u, z') f 2 (k, u,t) where 

h{k,u,t) = AeM-2(k 2 + u 2 )t] + 1 ~ eM ~ 2 2 ^ 2 + ^ (6) 

and 

f 2 (k,u,t)=Aex P [-(k 2 + u 2 )t}. (7) 
Therefore, the autocorrelation A(z,t) = J C s (k, z, z, t)d d ~ 1 k/(2ir) d ~ 1 is given by 

A(z,t) ~ r^ 2 ndueM-uH) {cSinU * + UC 2 ° SUz)2 . (8) 
Jo c 2 + u 

It is clear from Eq. (8) that in the limit z — > 00, we recover the bulk result: for large t, 
A(oo,t) ~ [£(£)] ~ d where £(t) ~ t 1 ^ 2 (Z = 2 within Gaussian theory) and hence Xb = d. 
On the other hand, for z = 0, we find that for large t, A(0, t) ~ [£{t)] ^^ d+2 ' > for c > and 
A(0, t) ~ [£(t)] ~ d for c = 0. Thus, we obtain the results that for the "special" transition 
(c = 0) A sp = d while for the "ordinary" transition (c > 0), X or = d + 2 within the Gaussian 
theory. It is interesting to note that while Xb satisfies the upper bound Xb < d conjectured 
by Fisher and Huse 0, clearly X or violates this upper bound. 

For d < 4, where the interaction term is no longer irrelevant, we evaluate the expo- 
nent A c in e = 4 — d expansion. The two-time correlator (unsymmetrized) C(x,x',t) = 
(4>{x', O)0(x, £)) in real space evolves as 

d t C(x,x\t) = [-ro + \7 2 }C(x,x',t) - -J2(M^ ,0)&(£, (9) 

n ij 

At the Wilson-Fisher fixed point, uq = 8ir 2 e/(n + 8) to leading order in e |IE |. This allows 
one to calculate the corrections to the two-point correlator perturbatively in Uq. To leading 
order in e, the term proportional to uq in Eq. (9) can be expressed, using Wick's theorem, 



4 



in terms of the mean field propagators as — uo(n + 2)Go(x, x, t)Co(x', x, t) where Go and Co 
denote the mean field equal-time and two-time propagators respectively. To leading order 
in e, one can replace Co in this term by C and then Eq.(9) becomes a linear evolution 
equation for the two-time correlator which is correct to 0(e). This evolution equation, for 
the symmetrized correlator C s (x,x',t), reads 

d t C s (x, t) = i[-2f + V £ 2 + V^ 2 - V(z, t) - V(z', t)}C s (x, t) (10) 

with r = r + Uo(n + 2)Go(oo, oo, oo) where Gq(x, x', t) denotes the mean-field equal-time 
propagator. The potential V(z, t) = uo(n + 2) J[Go(k, z, z, t) — Go(k, oo, oo, oo)]d 3 (fc)/(27r) 3 
captures the corrections due to fluctuations for d < 4 and can be calculated explic- 
itly. For example, at the "special" (c = 0) and the "ordinary" (c = oo) fixed 
points, we get, V sp (z,z,t) = — u (n + 2)[1 + [2t/z 2 ) exp(— z 2 /2t)}/32n 2 t and V or (z, z, t) = 
-ito (n + 2)[1 - (2t/z 2 ) exp(-z 2 /2t)]/32Tr 2 t in the scaling regime where z » A 1 , A being 
the upper cutoff. 

The Eq. (10) and the form of the potential V(z,t) suggest a late time scaling ansatz 
for the Fourier transform C s (k, z, z' ,t) rs t~ a exp(— k 2 t)f[z/\/t, z'/yt]. To determine the 
exponent a we first consider the bulk limit z — > oo, z' — > oo. Using V(oo,oo,t) = 
-it (n + 2)/32vr 2 t and u = 8vr 2 e/(n + 8) in Eq. (10) we get a = \ - (n + 2)e/4(n + 8) 
and f(x,y) ~ e~( x ~y) 2 as x,y — > oo. It follows immediately that the bulk autocorrelation 
Ab(t) ~ t~( d ~™+% z^ 2 . Since the dynamic exponent Z = 2 + 0(e 2 ), we recover the bulk 
result Xh = d — 2 r|l. For the surface autocorrelator, we need to know the small argument 

u n+8 2 ' ° 

behaviour of the scaling function f(x, y). For small x, y, f(x, y) ~ (xy) s (a + bx 2 + cy 2 + • ■ ■) 
where s(s — 1) = i^iff an d ± corresponds to "special" and to "ordinary" transitions re- 
spectively. Then the autocorrelator A(z,t) = J C s (k, z, z, ^d^k/ (27r) d - 1 ~ t -(d-i+2s+2a)/2 _ 
We choose the root of s to match the e — > limit and get X sp = d — -k f or "special" 
transition and A or = d + 2 — 2±| ^ for the "ordinary" one. 

We next calculate A c exactly in the large n limit. In this limit, C s (x,x',t) satisfies Eq. 
(10) exactly except that the potential V(z,t) is determined self-consistently from V(z,t) = 
uo(n + 2) J[G(k, z, z, t) — G(k, oo, oo, oo)]rf a!_1 (A;)/(27r) d_1 , where the equal-time propagator 
G(k, z, z', t) satisfies: 

d t G(k, z, z', t) = [-2k 2 + d 2 + d z , 2 - V(z, t) - V(z', t)]G(k, z, z', t) + 25(z - z'). (11) 

In analogy with epsilon expansion, we make the ansatz V(z,t) = ^ + ~^^ g[z/\/t\ where 
g(x) — > 1 as x — > and g(x) —>■ as x — > oo. The values of a and fi are determined 
respectively from the limits z — > oo (bulk dynamics) and t — > oo (static limit) and are 
already known to be a = (d — 4)/2 |13[, = (d— 5)/2 and /i or = (d — 3)/2 The 
full form of the scaling function g(x) is to be determined from a complicated self-consistent 
equation. However, for the purpose of determining the surface exponent A s , it is sufficient to 
know that g(0) = 1. We then proceed identically as in the case of e expansion by assuming 
a scaling ansatz for C s (k, z, z' ,t). We find a = (1 + a)/2 and s(s — 1) = fi 2 — 1/4. Using 
the fact that Z = 2 in the large n limit, we obtain X sp = (5d — 12)/2 and X or = (5d — 8)/2. 
Note that in the limit z — > oo we recover the bulk result Xb = (3d — 4)/2. These results are 
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consistent with those obtained from e expansion after taking n — > oo limit and also with the 
mean field results in d = 4. 

We have also carried out a direct numerical simulation of the two-dimensional Ising model 
with open boundary conditions at the bulk critical temperature using spin-flip Metropolis 
algorithm. We measure boundary spin correlations and compare them with the correspond- 
ing bulk measurements done with periodic boundary conditions. In the static limit, it is well 
known that the boundary spins order only due to the ordering of the bulk spins ( "ordinary" 
transition) |T4]-|T6||. Since at T c in the static limit, the boundary correlator falls off with 
distance r as 1/r (as opposed to the bulk decay r" 1 / 4 ) for large r |B|, it is natural to as- 
sume a dynamic scaling form for the equal-time boundary correlator G s (r, t) ~ r~ lr y(r /£(£)) 
at late times. It is believed that even in the presence of a boundary, there is still only a 
single time-dependent correlation length £(£) ~ t x l z with Z « 2.15 numerically pO| , pT| . The 
quantity Xs(t) = ((M s (t)) 2 ) (M s (t) being the total boundary magnetisation at time t) is the 
integral of the correlation function J s drG s (r,t) on the boundary. From the scaling form of 
G s {r,t), Xs(t) would then grow logarithmically with (and hence logarithmically with 
t). In contrary, the corresponding quantity in the bulk, Xb{t) = ((Mb(t)) 2 ) (Mb being the 
total bulk magnetisation), grows algebraically as [£(t)] 7 / 4 . In Fig. 1 we plot Xs{t) and Xb(t) 
against log(t) and, in the inset, show a log-log plot of Xb{t) versus t. Next, we compute the 
auto-correlation on the boundary of the open system and the periodic bulk system. The 
two results are contrasted in Fig. 2. The autocorrelation on the boundary decays much 
faster. This is expected for "ordinary" transition as evident from the large n and the e ex- 
pansion calculations. From the slope of the log-log plot in Fig. 2 we estimate the boundary 
exponent ratio A or / ' Z = 1.2 ± 0.1. The corresponding ratio in the bulk is estimated to be 



Xb/Z = .74 ± 0.02 in agreement with the previous simulation [12]. Using Z « 2.15, we 
estimate, X or rs 2.58 ± 0.1 to be contrasted with \ b m 1.59 ± 0.02. 

A two-dimensional system where these boundary auto-correlation exponents can be cal- 
culated exactly is the dynamics of the X-Y model following a quench from one temperature 
to another, both temperatures being below the Kosterlitz-Thouless temperature. In this 
case, the system is free of vortices and therefore, the dynamics of the phases is trivially that 
of damped independent spin waves. However, as shown by Rutenberg and Bray |J , the bulk 
auto-correlation decays algebraically with time and the corresponding exponent X b depends 
continuously on temperature. It is straightforward to extend their calculations to systems 
with edges and corners. For example, we have found that A corner = 2\ e d ge = 4A&. Details of 
these calculations will be published elsewhere W2 . 



It is clear from the results presented in this paper that for quench to T = T c , the 
surface autocorrelation exponents are quite different from their bulk counterparts. It is 
an interesting question to ask whether the same happens for quench to the ordered phase 
(T < T c ). From our preliminary T = simulations of the 2d Ising model, A s does not seem 
to be different from A& at T — 0. This seems to be the case also for the exact T = Glauber 
dynamics of the semi-infinite Ising chain. Interpretation of these results and further studies 
are relegated to a future publication 
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FIGURES 



FIG. 1. Xs(t) is plotted against log(i). The logarithmic dependence is pretty evident. The 
inset shows a log-log plot of Xb(t) versus t, which is consistent with a power law. 

FIG. 2. The autocorrelation on the boundary, A s (t) versus time t in a log-log plot. From the 
slope we estimate the boundary exponent ratio A or /i? = 1.2 ± 0.1. The corresponding ratio in the 
bulk is estimated to be \hj Z = .74 ± 0.02. Bulk data is shown in the inset. 
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